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SUMMARY . 



This report has been prepared to provide a' practical 
method for determining the chordwise distribution of the 
rate of heat transfer from the surface of a wing or body 
of revolution to air. The method is limited in use to 
the determination of heat transfer from the forward sec- 
tion of such bodies when the flow is laminar. A compar- 
ison of . the calculated average heat-transfer coefficient 
.for the nose section of the wing of a Lockheed 12-A air- 
plane with that experimentally determined shows a satis- 
factory agreement. A sample calculation is appended,. 



INTRODUCTION 



With the advent of thermal ice-prevention equipment 
of aircraft, the problem of determining the heat which 
may be transferred from the surface of aerodynamic bodies 
has become of considerable interest. Experimental inves- 
tigations (reference l) have indicated that adequate heat- 
ing of the forward 10 to 20. percent of an airfoil will 
prevent ice formation on the entire airfoil. Moreover, 
experience has shown that in the cruising flight condition 
in the absence of • icing- conditions if sufficient heat can 
be supplied' to this section of the airfoil to raise its 
surface temperature from- 70° t © 100° F above the tempera- 
ture of the ambient air, ice will not' collect on the air-, 
foil. 

To determine whether this surface temperature rise 
can be obtained in the design of a particular installa- 
tion, it is necessary to determine the rate of heat trans- 
fer from the airfoil surface to'the atmosphere. In the 
usual crui s ing-f 1 ight condition the nose section of an 
airfoil experiences laminar flow. In the present report 
a method for calculating the rate of heaV transfer from 
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an airfoil section subject to laminar flow is developed. 
The method is made general to include laminar flows oc- 
curring over surfaces subjected to positive or negative 
pressure gradients and may be applied to two-dimensional 
bodies as well as bodies of revolution. 



The following development of an expression for the 
heat transfer between solids and fluids has been based 
upon Reynolds analogy. This analogy was presented by 
Reynolds in an early paper (reference 2) in which he- sug- 
gested that, in a fluid, momentum and heat are transferred 
in the same way, and concluded that in geometrically sim- 
ilar systems a simple proportionality relation exists be- 
tween heat transfer and fluid friction. 

The analogy may be applied in the case of th+e heat 
transfer between solids and fluids when the fluid flow is 
laminar, .provided the Prandtl number for the fluid is 
unity. The Prandtl number is the dimensionless parameter 



c specific heat of the fluid at constant pressure 



k thermal conductivity of the fluid 

Tor air the value of the Prandtl number is 0.73 
rather than unity but, as is discussed later- in this 
report, experimental investigations have shown that heat 
transfer is only slightly affected by variation of the 
Prandtl number so that for most practical purposes the 
Reynolds analogy may be applied. The analogy leads to 
the equation 



THEORY 



Pr = 




where 



absolute viscosity of the fluid 



h = 2 c f pc p V 1 



(1) 



where 



h heat transfer coefficient defined as the heat trans- 
ferred per unit time from a unit surface area for 
1° difference of temperature "between the surface 
and the fluid "outside 11 the boundary layer 

p fluid density 

V 1 stream velocity t just outside the boundary layer 

Cf surface friction coefficient defined as the frictional 
drag per unit surface area in terms of the local 
fluid dynamic pressure, q a 

The surface friction' coefficient may he expressed as 



where T is the frictional force per unit area. 

For laminar flow the frictional force T is the 
product of the absolute viscosity and the velccity gradient 
in the fluid boundary layer at the solid surface. For any 
•given b oundary-^layer * vel ocity profile, this velocity gra- 
dient is directly proportional to the velocity Y 1 and 

inversely proportional to the boundary-layer thickness. 
It follows that . 

where 5 is the b oundary-lay er thickness and A is a 
constant dependent on the shape, of the boundary-layer 
velocity profile and on the definition of. 5. It has be- 
come customary (reference 3) to define 6 for lamiriar- 
boundary layers as th? distance from the solid surface to 
a point in the boundary layer where the dynamic pressure 
is one-half its value outside the boundary layer. 

Inserting the value, of Of from equation (2b) into 
equation (l) the heat-transfer coefficient becomes 
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Rearranging and dividing "both sides by the thermal 
conductivity of the fluid k this "becomes 



X Pr = X (3). 

since the Prandtl number must he assumed to he unity by 
the analogy. The parameter (h5)/k is hondimensional 
and might be properly termed "the boundary-layer Nusselt 
numb er • H 

The appearance of Pr in equation (o) might suggest 
that its value for air (0.73) should be substituted to ob- 
tain a better approximation of the 'value of the boundary- 
layer ITusselt number. However, the meager exper iment-al 
evidence available (reference '4 , p. 249) indicates that 
for- laminar flows, the Nusselt number is proportional to' 

(Pr) 4 so that using. unity for. the value of Pr (equiv^' 
alent to (Pr).) would appear to be preferable. 

In order to determine h' from equation (?) it is 
necessary to know the value of A which is a'function 
only of the shape of the boundary-layer velocity profile. 
Numerous experiments with airfoils have shown that the 
velocity profile of the laminar boundary layer in the 
presence of favorable pressure gradient (i.e., where 
the surface pressure gradient is negative proceeding in. 
the downstream direction) is closely approximated .by the 
Blasius distribution for the flat plate. For the . B.l.as.ius 
type, velocity profile the value of X is 0.765. 

At points on the surface of a body downstream of the 
minimum-pressure point a laminar boundary layer exhibits 
a tend en cy to separate. Since the velocity gradient at 
the surface decreases as separation devel ops the value of 
X must diminish' until at the separation point its .value 
is zero. It is considered that reducing A linearly 
from the minimum pressure point to the separation point 
will satisfactorily approximate the actual pase, particu- 
larly since this region is of little importance from the 
viewpoint of heat transfer. The location of the laminar 
separation point may be calculated by the method of ref- 
erence 5 . 
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With the known variation of X along the surface, 
the heat- trans fer coefficient may be calculated when the 
boundary-layer thickness has teen determined. The value 
of 6 may be found for the region from the stagnation 
point to the minimum pressure point" by the method of ref- 
erence 3. Rearranged in a form more convenient for calcu- 
lation, for airfoils the equation is 



8° = 



5.. 3c' 



Sj/c 

/ 

•o vv o 



6.17 



L 



where 



(4) 



s 

«i 

▼■o 
V 

R, 



airfoil chord length 

distance along the surface from the stagnation point 

distance to- the point for which the "boundary- layer 
thickness is "being computed 

free- stream v eloci-ty ■■ 

velocity outside the ■boundary layer at s/c . _ 

velocity outside the "boundary layer at Sj/c 

cV rt 

"body Reynolds number = 



V 



For "bodies of revolution the corresponding equation is 




6 



where, in addition. 
L b ody length 

r radius of revolution at s 
r a radius of revolution at Bj 



body Reynolds number = 



The value of 6 is indeterminate at the stagnation 
point. It" is not necessarily zero as may he -shov/n by the 
following approximate analysis. The initial flow near the 
stagnation point on an airfoil will he approximately that 
over a circular cylinder having a radius equal to the 
radius of curvature at stagnation F, On the surface of a 
circular cylinder the velocity at 6, the angular coor- 
dinate measured from stagnation, is 



JL 
V 



= 2 sin 0 



For small values of 9 this becomes 



= 28 = 



r 



(r/c) Vc 



Using this value in equation (4), gives 

s a /c 



6" = 



2. 65c 1 - 



(0 



9.17 



3.17 



so that for small- values. 'we' obtain the approximate value. 



' 6 



stag 



0, 28c 3 F 
R_ % c 



A better approximation is to ddnsider the nose of the air- 
foil ellpitical in form. By such an approximation, if 
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the stagnation point is located at the leading edge, the 
calculated value of 5 s tag lies between the value obtained 
for the cylindrical nose and half this value, depending on 
the thickness ratio. Neither the order of approximation, 
nor the importance of the value of S stag warrants much 

refinement of this calculation* It is recommended that 
the .approximation 

; 6 3 «!LE - ( 6 ) 

•stag 5R C c 

he used for two-rdimens i onal bodies,, where F is the radius 
of curvature at the 3 1 agnat i on . point \-/hich is not neces- 
sarily the leading-edge radius. ' 

For bodies of revolution, by a similar- treatment it 
may be shown that equation (6) is, with; ; R c replaced by* 

R L , a fair approximation if r is the least radius of 

curvature' at the stagnation point. 

As "has been noted, the equations (4} and (5) were 
devised on the assumption that the vel oci ty. prof i le of 
the boundary layer was of the Blasius type. In the region 
behind the minimum' pressure .point the velocity profile 
changes to a separated profile so that ' if these equations' 
were used to determine '6 in this regi.on, it would be 
expected' that the calculated, would exceed -the actual 
boundary-layer thickness. A comparison of the. experi- 
mentally determined boundary-layer thickness for the NACA 
0012 airfoil (reference 6) with that calculated, by this., 
equation and by the more 'exact but laborious method of . . : 
reference (7) shows' (see fig. l) that this equation 
* appears to yield satisfactory results .over tire separating 
region up to' the separation point. 

In the preceding theory tie effect of the temperature 
distribution within the 'heated boundary layer on the. 
boundary-layer thickness and velocity profile has been 
neglected. This effect is known 'to be small .for the tern-'- 
perature differences necessary for thermal ice-prevention 
and need not be c ons idered .f or such practical applications. 
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DISCUSSION AND CONCLUSIONS 



Recently an experimental investigation of a heat de- 
icer installation on a Lockheed 12-A airplane wing in 
flight was conducted. In the course of this investiga-.. 
tion'the average heat-transfer coefficient from the for- 
ward 10 percent of the airfoil surface to the air was 
determined. 

For purposes of comparison the chordwise distribution 
of the heat-transfer coefficient was calculated for the 
Lockheed wing at three spanwise stations and the average- 
heat-transfer coefficient determined. The calculated 
value of 13 Btu per hour per square foot per degrees Fahr- 
enheit agrees well with the value of 11 Btu per hour per 
square foot per degrees Fahrenheit as determined from 
flight tests. 

■. In the past, for lack of appropriate experimental 
data, heat-transfer coefficients have been calculated on 
the assumption that the experimental data of reference 1 
for the Clark Y airfoil at a given angle of att'ack may 
be considered to apply to airfoils of other shapes pro- 
vided the results are corrected for the effect of scale 
found to' apply . over thelimited range of test Reynolds- 
numbers for the Clark Y airfoil. On this basis, the 
heat-transfer coefficient for the Lockheed 12-A airplane 
wing would be 18 Btuper hour per square foot per degree 
Fahrenheit. The agreement ; such as it is, must be. con- 
sidered fortuitously close. The assumption that, .at' the 
same angle of attack and Reynolds number, the heat- 
transfer coefficient for the Lockheed wing (NACA 230^ 
se'ries sections) and the Clark Y will b.e identical is 
clearly unjustified.- Moreover, there is no. basis far the 
assumption that the heat-transfer characteristics at high 
Reynolds" numbers may be determined by extrapolation of 
low Reynolds number test results using an extrapolation 
formula of the kind 

> = » ! < 

where' n and a are constants. 

Since the heat-transfer coefficient is determined by 
the boundary-layer thickness which in turn is dependent 
on the pressure distribution, it is essential that the 



pressure distribution used in" th-e determination of the 
heat-transfer coefficient by the method of this report, 
must he that corresponding to the proper wing Reynolds 
number. In the event that the required experimental ■ 
pressure distribution is not available,., the use of the 
methods of references 8 and 9 combined or the more labo- 
rious but more exact method of reference 10 as modified, 
by thie method' of reference 11 is recommended, Without 
this modification, the method of reference 10 isunsat- • 
isfactory at high angles of attack as is indicated in 
refer ence 11 • 

In regard- to the application of the calculated wing 
surface-to-air heat-transfer coefficient given by equa- 
tion (3) to the determination of the rate of heat trans- 
fer from the airfoil surface, it should be noted that 
the temperature difference used in this calculation 
should be the surface t emperature. minus the air temper-, 
ature as it is ' increased as a result of fluid friction. 
The air temperature rise due to fluid friction is ap- 




proximately l* 7 V77rkj in degrees Fahrenheit, where 7 



0 



is the airplane airspeed in miles per hour so that the 
rate of heat transfer is 

where 

Q, rate of heat transfer 

S heated area 

T s surface temperature 

T 0 ambient air temperature 

and 'all variables are in consistent units. The correc- 
tion due 1 0 ' aerodynamic Heating is only important at 
very, high speeds and so normally may be neglected. 
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APPENDIX 



Sample calculation: 

In order to illustrate the method of determining the 
chordwise distribution of the heat-transfer coefficient 
developed in this report., calculations were made for the 
wing of the Lockheed 12-A airplane used in the experimen- 
tal flight investigation of heat de-icing. The wing sec-, 
tion used for this illustration was taken at a spanwise 
station 123 inches from the center line of the fuselage, 

For this calculation the following were assumed: 

Air temperature, °F . . . . 30 

Altitude, feet 8125 

True airspeed, miles per hour 173 

The wing-lift coefficient was 0.424. At spanwise* 
station 123 the chord was 7,78 feet and the airfoil thick- 
ness was 14 1/ 2 percent of the chord. Reynolds number, 
based upon the chord length, was 

R c s = 10,330,000 

where 

c chord (7.78 ft ) 

V velocity of ambient stream (253 ft/sec) 

V kinematic viscosity (0.0001906 ft a /sec) 

It was first necessary to calculate the chordwise 
pressure distribution because no data* could be found 
which were applicable to this case. The section-lift 
coefficient was computed by the method given in- reference 
12. The chordwise normal-force distribution was then 
found, using the method given in reference 8. To find 
the chordwise .pressure distribution on the upper and 
lower surfaces of thr- airfoil, and from this the corre- 
sponding velocity distribution., the method of reference 3 
was used. These calculated values have been tab\ilated in 

the form (l-P) and r^, where ' V 0 ' is the'velocity in 

v o 
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ambient stream. (See table I.) A plot of the pressure 
distribution was made to determine the stagnation point 
(fig. 2) the position of the laminar-separation point on 
the .upper surface was found by the method of reference 5, 
(See fig. 3 . ) 

The heat-transfer coefficient h was calculated by 
equation (3) of this report 

h = *A 

6 

6, the laminar boundary-layer thickness was computed by 
equation (4) of this report 



5.3c' 



a /c 



e . i 7 



Substituting value of c = 7,78 feet and R Q 
10,330,000, this becomes 



= 0.000031 



e . i 7 






was evaluated graphically, and is given in table I. 
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Equation (6) was used to find 6 near the stagnation 
point. In the calculation of h the value of X used 
was 0.7.65, except between the minimum pressure point and 
the laminar-separation point on the upper surface of the 
airfoil. In this region X was reduced linearly from . 
0.765 at the minimum pressure point to 0 at the laminar- 
separation point. This variation of X has, "been dis- 
cussed in the section of the report giving the theoretical 
development of h. 

For these calculations 

1.24 x 10~ 2 X 

h = — — i 



Substituting for k = thermal conductivity of the air at 
30° F = 1.24 x 10~ 2 Btu per hour per spare foot per 
degree Fahrenheit per foot. 

Values of X, 5, and h are given in table I. 

The chordwise distribution of 5 and h have been 
plotted in figure 3 • 
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Figure 1.- The calculated and measured "boundary layer thick- 
ness on the surface of an NACA 0012 airfoil at 

zero lift, 



BACA 



FigQ. 3,3 



3.0 



1.6 



1.2 



U-P) 



0.8 



0.4 



A 



Stagnation point 



P - P, 



■°c pressure coeficient 



Po f 



static pressure on 
in 



q = dynamic 
-(1-Py) U refers to upper surface 



-H 



M1-P L > 



airfoil surface 
ambient stream 



Loner 



8 12 16 20 24 

Distance from nose of airfoil, percent chord, JL xlOO 



32 



Figure 2.- The chordwise pressure distribution at ning station 133 of the Lockheed 
12A airplane, at wing C L = 0.424 and Rc = 10,330,000. 
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Figure 3.* Velocity boundary layer thickness and heat tranofer coefficient distribution 
along the surface at wing station 133 of the Lockheed 13 A airplane at 
wing C L « 0.434, 1^ « 10,330,000. 



